Multiple Sources toward the High-mass Young Star S140 IRS 1 



Miguel A. Trinidad^ Jose M. Torrelles^, Luis F. Rodriguez'^, Salvador Curiel^ 

trinidadOastro . ugto . mx 

ABSTRACT 

S140 IRSl is a remarkable source where the radio source at the center of the 
main bipolar molecular outflow in the region is elongated perpendicular to the 
axis of the outflow, an orientation opposite to that expected if the radio source 
is a thermal jet exciting the outflow. We present results of 1.3 cm continuum 
and H2O maser emission observations made with the VLA in its A configuration 
toward this region. In addition, we also present results of continuum observations 
at 7 mm and re-analyse observations at 2, 3.5 and 6 cm (previously published). 
IRS lA is detected at all wavelengths, showing an elongated structure. Three 
water maser spots are detected along the major axis of the radio source IRS 
lA. We have also detected a new continuum source at 3.5 cm (IRS IC) located 
~ (y.'6 northeast of IRS lA. The presence of these two YSOs (IRS lA and IC) 
could explain the existence of the two bipolar molecular outflows observed in the 
region. In addition, we have also detected three continuum clumps (IRS IB, ID 
and IE) located along the major axis of IRS lA. We discuss two possible models 
to explain the nature of IRS lA: a thermal jet and an equatorial wind. 

Subject headings: YSOs — ISM: individual (S140 IRS 1) — ISM: jets and outflows 
— masers — stars: formation 
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1. Introduction 

The process of low-mass star formation is reasonably well understood. The accepted 
model (e.g., Shu et al. 1987) requires a disk-YSO-outflow system ( YSO: young stellar object), 
and it has been supported by theoretical and observational results (e.g. Evans, 1999 and 
references therein). In addition, bipolar molecular outflows associated with low-mass stars 
are believed to be driven by their jets, that can be observed at sub-arcsec scales and this 
constitute the best evidence of coUimation at the smallest scales now known (e.g. Anglada 
1996). However, the process of high-mass star formation is not well understood yet. In 
fact, although molecular outflows seem to be also common among high-mass stars (Gomez 
et al. 1999, Zhang et al. 2001, Ridge & Moore, 2001; Beuther et al. 2002, Shepherd 2005), 
there is still a deficit in the detection of circumstellar disks and jets associated with massive 
YSOs (e.g., Cepheus A HW2: Rodriguez et al. 1994; Patel et al. 2005; Curiel et al. 2006; 
IRAS 20126+4104: Cesaroni et al. 1999; Trinidad et al. 2005; Sridharan et al. 2005; AFGL 
490: Schreyer et al. 2006; G24.78+0.08: Beltran et al. 2006). In spite of these studies, it 
is not clear whether molecular outflows in massive YSOs are also driven by coUimated jets 
as low-mass stars do. Therefore, studies of individual high-mass YSOs are very important 
to address these issues, but their identification is more difficult than in the case of low-mass 
objects. Firstly, the spatial density of massive objects is smaller and we need to go to greater 
distances to increase the known sample and build statistical significance, and second, they 
are formed in dense clusters and sub-arcsecond angular resolution observations are required 
to isolate them. 

The star-forming region, S140 IRS, is located in the molecular cloud L1204 at a distance 
of 910 pc (Crampton & Fisher, 1974). In this region, there is a group of at least three ZAMS 
B stars (IRS 1, IRS 2, IRS 3; Beichman et al. 1979), with a total luminosity of ~ 3 x 10^ 
(Lester et al. 1986), with IRS 1 being the brightest one. In addition, a CO bipolar outflow 
extended along the northwest-southeast direction has been observed in the region by Hayashi 
et al. (1987). Minchin et al. (1993) found that the blue and redshifted lobes of the CO 
bipolar outflow have position angles of ~ 160° and ~ 340°, respectively. On the other hand, 
K band (2.0-2.3 /im) and H2 observations have revealed two bipolar outflows in the region 
(Preibisch & Smith 2002 and Weigelt et al. 2002), one of them with an orientation similar 
to the CO outflow (160/340°) and the other one in the 20/200° direction. Both bipolar 
outflows seem to be centered on IRS 1. The region S140 IRS has also been studied at the 
infrared and optical bands (Eiroa et al. 1993), NH3 lines (Verdes- Montenegro et al. 1989), 
and radio continuum emission (Schwartz 1989, Evans et al. 1989). More recently, Hoare 
(2006) has observed IRS 1 at 6 cm during three epochs with MERLIN, and shows, through 
a very detailed analysis that it is highly elongated in the northeast-southwest direction, 
and also proposes that the radio continuum traces an ionized equatorial wind driven by 
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radiation pressure from the central star and oriented in the northeast-southwest direction, 
perpendicular to the CO bipolar outflow. In addition, water maser emission has also been 
detected toward the S140 IRS region (e.g. Tofani et al. 1995, Lekht et al. 1993, Lekht & 
Sorochenko 2001, Trinidad et al. 2003). 

In this paper, we report and discuss new high angular resolution observations of 1.3 cm 
continuum and water maser emission toward the S140 IRS region carried out with the Very 
Large Array (VLA) in the A configuration. In order to present a full study of IRS 1, we have 
also analyzed 7 mm continuum observations (VLA data archive; see also Gibb & Hoare, 2007) 
and re-analyzed VLA centimeter continuum observations (previously published by Schwartz 
1989 and Tofani et al. 1995). We describe the observations in §2 and present the results in 
§3. In Section §4, we discuss the nature of the multiple continuum sources that we detect in 
the region, while our main conclusions are summarized in §5. 



2. Observations 

The observations toward the S140 IRS region were made with the VLA of the Na- 
tional Radio Astronomy Observatory (NRAO)Q in its A configuration on 1999 June 29. We 
observed simultaneously 1.3 cm continuum and water maser emission. We used two differ- 
ent bandwidths, one of 25 MHz with seven channels for the continuum and another one of 
3.125 MHz with 63 channels for the maser line emission. Both the right and left circular 
polarizations were sampled in the above two different bandwidths, which were averaged in 
order to improve the sensitivity. The broad bandwidth for the continuum observations was 
centered at 22285.080 MHz, while the narrow bandwidth for the line observations was cen- 
tered at the frequency of the H2O Gig 623 maser line (rest frequency 22235.080 MHz) with 
Vl5/j = —9.0 km s~^, covering from —28 to 10 km s~^ in velocity. The absolute amplitude 
calibrator was 3C 286 with an adopted flux density of 2.51 Jy, while the phase calibrator 
was B2021+614 with a bootstrapped flux density of 2.14 ± 0.05 Jy. The water maser line 
and continuum data were reduced and calibrated using the standard techniques using the 
NRAO AIPS software package. After the first calibration, we searched the narrow bandwidth 
data for the spectral channel with the strongest water maser emission and then its signal 
was self-calibrated in phase and amplitude. We then cross-calibrated the data applying the 
phase and amplitude corrections to both the narrow and the broad bandwidths. In this way, 
the atmospheric and instrumental errors at high frequencies on long baselines were removed 
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(see for details Reid & Menten 1990, Torrelles et al. 1996) and the signal-to-noise ratio was 
improved. 

We also used the data set at 7 mm from the VLA data archive (see also Gibb & 
Hoare, 2007). In addition, we re-analyzed continuum observations at 2, 3.5 and 6 cm, which 
have been previously published (see Tabled]). All observations were carried out in the A 
configuration. We recalibrated these data sets and made new contour maps using the current 
procedures of AIPS. We used the physical parameters derived from these new contour maps 
for the discussion in this paper. Contour maps at 2 and 6 cm and their measured physical 
parameters are consistent with those published before by Schwartz (1989). However, the 
flux density estimated at 3.5 cm is different (about 50% lower) to that reported by Tofani 
et al. (1995), probably due to a typo. 

3. Observational Results 

Figure [1] shows the contour maps of S140 IRS 1 at centimeter wavelengths (1.3, 2.0, 
3.5 and 6.0 cm). Given that other condensations are also detected around S140 IRS 1, 
hereafter, we will refer to S140 IRS 1 as S140 IRS lA or only IRS lA for clarity. This source 
appears spatially resolved at all wavelengths. The physical parameters of the radio source 
IRS lA (position, flux density, deconvolved angular size and position angle) were obtained 
from elliptical Gaussian fits using the AIPS task IMFIT and are given in Table [2l Contour 
maps at 1.3 and 2 cm have similar angular resolution (~ O'.'l) and were made with natural- 
and uniform-weighting, respectively. On the other hand, contour maps at 3.5 and 6 cm have 
an angular resolution of ~ 0'.'3 (Figure [T]). For the map at 3.6 cm a Gaussian taper of 850 kA 
was used, while at 6 cm a ROBUST = -1 parameter (Briggs 1995) of the AIPS task IMAGR 
was used. 

In all contour maps (Figured]), IRS lA shows a general elongated morphology in the 
northeast-southwest direction, similar to the orientation of the bipolar outflow observed in 
Ha toward IRS lA in the 207200° direction (Preibisch & Smith 2002). In addition, our 
contour map at 1.3 cm with angular resolution of ~ 0'.'08 (ROBUST = parameter was 
used in order to optimize the compromise between the signal-to-noise ratio and the angular 
resolution) shows a second peak, at a level of 4a, located about 0'.'13 to the southwest from 
IRS lA (Figure[2]). We will refer to this peak as IRS IB. Furthermore, three other continuum 
peaks (that we will refer as IRS IC, IRS ID and IRS IE), are also observed at 3.5 cm with 
an angular resolution of 0'.'25 (Figure [2]). IRS IC is located ~ 0'.'6 to the northeast of IRS lA 
(Figure [2]) and is also detected at 2 cm with an angular resolution of ~ 0'.'15. IRS ID and 
IRS IE are located ~ 0'.'35 to the northeast and southwest of IRS lA, respectively. Their 
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main observed parameters are given in Table [2l 

On the other hand, IRS lA is detected at 7 mm (Figure [31 flux density and deconvolved 
size are given in Table [2]) while this is not the case for IRS IB, IC, ID or IE. Although 
the 7 mm continuum emission is also elongated in the northeast-southwest direction as the 
centimeter emission (1.3, 2, 3.5 and 6 cm), they do not have, within the error, the same 
position angle (the millimeter emission has a position angle of ~ 61°, while the centimeter 
emission has a position angle of ~ 44°). These results are consistent with those found very 
recently by Gibb and Hoare (2007). 

The water maser emission toward S140 IRS was in a period of minimum activity during 
our VLA observations and we only detect three water maser features spatially associated 
with IRS lA (Figured]). One of these maser features, with redshifted velocity with respect 
to the molecular cloud velocity (~ —6.5 km s~^, Zhou et al. 1993), is located close (O'.'l) to 
the IRS lA continuum emission peak at 1.3 cm, while the other two maser features, with 
blueshifted velocity with respect to the molecular cloud velocity, are located ~ 0'.'5 and ~ 0'.'7 
to the southwest of IRS lA, respectively (see Fig. [T]and Table [3D • Figure 9 from Trinidad 
et al. (2003) shows that the water maser emission toward S140 IRS was stronger (up to 46 
Jy with at least seven features) a month before the VLA observations presented here. The 
position, velocity, and flux density of the three detected VLA water maser features are given 
in Table [31 

4. Discussion 

The elongated morphology of the IRS lA is evident in the contour maps at all centimeter 
wavelengths (1.3, 2, 3.5, and 6 cm), suggesting a jet-like nature (see Figure [T|), similar to 
other possible radio jets associated with massive YSOs (e.g., Claussen et al. 1994, Rodriguez 
et al. 1994, Torrelles et al. 1997, Hofner et al. 1999, Gibb et al. 2003, Trinidad et al. 2003, 
2005, Curiel et al. 2006; see Hoare 2006 and references therein). However, Hoare (2006), 
based on very detailed 6 cm continuum multiepoch observations with MERLIN, has proposed 
that the continuum emission is produced by an equatorial wind. We discuss both of these 
possibilities below. 

4.1. IRS lA: A Thermal Radio Jet? 

In order to investigate whether or not IRS lA is a thermal jet, we have calculated, 
following the formalism of Reynolds (1986), the dependence of the flux density and the 
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deconvolved angular size of the major axis of the jet with the frequency. For the case of a 
thermal jet with constant velocity, temperature, and ionization fraction, the dependence of 
the flux density with the frequency is Sy oc z/" (a = 1.3 — 0.7/e; equation 14 of Reynolds 
1986), while that for the deconvolved angular size of the major axis is Omaj oc z/^ (/5 = 
—0.7/e). Under this formalism, e is the power-law index that describes the dependence of 
the jet half- width, w (perpendicular to the major axis of the jet), as a function of the distance 
from the origin. For a conical jet, e = 1, then, a = 0.6 and (3 = —0.7. 

From the multiepoch observations at 6 cm (Hoare 2006), it is known that the flux 
density of IRS lA is not variable with time. Assuming that the flux density of IRS lA has 
not changed significantly with time also at other wavelengths (the observations were made 
with a time span of about 12 years; see Tabled]), we roughly estimate a spectral index of 
a = 0.5 ± 0.1 by using the flux densities at 0.7, 1.3, 2, 3.6, and 6 cm (see Figure HI). In 
addition, we also estimate the dependence of the angular size of the major axis of IRS lA 
with frequency, where the parameter is /3 = 0.98 ± 0.06 (see Figure H]). Both indices, a and 
/3, seem to be consistent with free-free emission in the wavelength range 0.7-6 cm, produced 
either by a conical coUimated ionized jet (e = 2/3, "standard" collimated jet) or by a conical 
ionized jet (e = 1, "standard" spherical jet) (see Table 1 of Reynolds 1986). On the other 
hand, using the flux density of IRS lA at 1.3 cm (this work) and that at 6 cm obtained 
with MERLIN observations (Hoare 2006), a spectral index of about 0.6 ±0.1 is obtained, 
which is consistent with that estimated using only VLA observations. In this way, we also 
note that the flux density and angular diameter of IRS lA as measured with the VLA, with 
an angular resolution of about 0'.'35, and those measured with MERLIN (Hoare 2006), with 
an angular resolution of about O'.'ll, are consistent within 25-30%. This result shows that 
changes of the angular diameter of the source is not due to an effect of angular resolution; 
that is, the larger angular size measured at low frequencies is not an effect of lower angular 
resolution than at higher frequencies. 

Under the assumption that IRS lA is a thermal jet, we can also estimate the ionized 
mass-loss (M) and momentum rate (P) deposited by IRS lA into the ambient medium. 
Following Reynolds (1986) and using equation 3 given by Beltran et al. (2001), for a pure 
hydrogen jet with constant velocity and ionization fraction, we have 
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where a is the spectral index, 5*;^ is the observed flux density at frequency z/, V is the terminal 
velocity of the jet, is the turnover frequency, 9o is the opening angle, i is the jet axis 
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inclination, d is the distance to the source and T is the electron temperature. The opening 
angle is approximately estimated using 9o = 2ta.n'^{9„iin/(^maj)- 

Then, assuming a inclination angle of the jet axis relative the line of sight, i = 60°, a 
distance of 910 pc, an electron temperature of 10'^ K, a lower limit for z/^ (= 22.2 GHz), 
and a terminal velocity of ~ 500 km s~^, we determine the mass-loss (M) and momentum 
rate (P) deposited by IRS lA into the ambient medium as ~ 9 x 10^^ M0 yr^^ and ~ 
4.5 X 10~^ Mq yr~^ km s~^, respectively. These values are similar to those obtained for 
the thermal jets associated with other massive YSOs (e.g Cep A W2 and HH 80). On the 
other hand, using the correlation (momentum rate versus radio continuum luminosity) given 
by Anglada (1996) to estimate the momentum rate, we find a value about two order of 
magnitude larger than that estimated using the formalism of Reynolds. This result suggests 
that Anglada's correlation is valid only for low luminosity objects and cannot be extrapolated 
to high luminosity sources. 

We have also detected four condensations almost aligned along the major axis of the 
IRS lA (see Figure[2]). One continuum peak (IRS IB) is detected at 1.3 cm, while the other 
three continuum peaks (IRS IC, ID, IE) are detected at 3.5 cm. As IRS IB is detected only 
at 1.3 cm (although IRS IB and IE are located to the southwest of IRS lA, both continuum 
peaks are not spatially coincident), we are not able to estimate its spectral index, which could 
have helped to investigate whether it is a condensation ejected by IRS lA or, alternatively, an 
independent source. The continuum peaks IRS ID and IE are almost symmetrically located 
with respect to IRS lA (see Figure El) and both condensations have similar peak flux densities, 
which could suggest that they are condensations ejected by IRS lA, as could be expected for 
a thermal jet. Multiepoch observations (e.g at 3.5 cm) to measure possible proper motions 
of IRS ID and IE could test if these are condensations ejected by the jet. On the other 
hand, the condensation IRS IC does not have a counterpart on the symmetrically opposite 
side as in the case of the pair formed by ID and IE. In addition, IRS IC appears misaligned 
with respect to the ID-IRS lA-lE system (see Figure [2]). Then, if IRS IC was also ejected 
by the jet IRS lA, this could indicate that IRS lA is precessing. Alternatively, IRS IC 
could be an independent source. In this way, IRS IC seems to coincide spatially, within 
the positional error, with the clump NE-N detected by Hoare (2006) at 6 cm wavelength. 
Using the peak flux densities at 3.5 (this work) and 6 cm (Hoare 2006), we have roughly 
estimated a spectral index for IRS IC of ~ 0.4, which suggest that the continuum emission 
of IRS IC is being produced by partially-thick thermal free-free emission. In addition, Hoare 
(2006) has reported possible proper motion of the NE-N clump with a tangential velocity of 
120 km s~^ moving toward the east. However, due to this high value of the velocity, Hoare 
(2006) opens the possibility that the proper motion is not due to real physical motion but 
due to a travelling illuminating wave. 
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The overall radio continuum properties of IRS lA indicate that it could be a thermal jet 
associated with a high-mass YSO. In addition, the position angle of the elongated structure 
at cm wavelengths (~ 44°; Table [2]) suggests that IRS lA could be the driving source of 
the bipolar molecular outflow seen in H2 in the direction 20°/200°. Under this scenario and 
assuming that the condensation IRS IC is an independent source, we speculate that this 
newly detected source (IRS IC) could be the driving source of the other bipolar molecular 
outflow observed in the direction 160°/340°. If this is the case, two close, independent and 
almost perpendicular jets could be the driving sources of the two bipolar outflows observed 
in the region (see Figure [2]). Perpendicular thermal jets have been observed in low-mass star 
forming regions (e.g. HH 111, Reipurth et al. 1999). 

In addition, all water maser features detected toward S140 IRS are associated with 
IRS lA, assuming that IRS IE is not a stellar object but ejected from IRS lA (see Figured] 
and [2]). The estimated mass for IRS lA, assuming that the masers are bound gravitationally, 
is > 60 M0, which is greater than the mass estimated from its luminosity (~ 10 M©). This 
result indicates that the water maser features are not associated with a circumstellar disk, 
instead, the water masers are tracing unbound motions, which is consistent with IRS lA 
being a thermal jet. Further observations with high angular resolution and proper motion 
studies will be needed to study in detail the three-dimensional (spatio-kinematical) distri- 
bution of the water masers in the region, which will help to understand the nature of the 
IRS lA and possibly discriminate between the thermal jet and equatorial wind models (see 
below) . 



4.2. IRS lA: An Equatorial Wind? 

Using high angular resolution multiepoch observations at 6 cm toward IRS lA, Hoare 
(2006) proposes that IRS lA is an equatorial wind driven by radiation pressure from a central 
star and inner disk acting on the gas in the surface layers of the disk, and perpendicular to 
the CO bipolar outflow observed in the direction 160°/340°. This scenario is supported by 
two facts. First, the location and extension perpendicular to the major axis of IRS lA of 
a small scale monopolar near-IR reflection nebula at the base of the blueshifted lobe of the 
molecular outflow, and second, the small outward proper motion of a continuum clump in 
the direction 160°/340° as would be expected for a jet. 

The results of Hoare are opposed to the picture of the multiple jets described above 
for IRS lA. However, our observational results could be also consistent with IRS lA being 
an equatorial wind. As mentioned above (§ 4.1), there are two condensations (IRS ID and 
IE) located symmetrically on opposite sides of the central source, IRS lA, and under the 
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equatorial wind model, IRS ID and IE could be tracing the outer border of the wind within 
a disk with diameter of ~ 650 AU. 

Unfortunately, with the present data we cannot discriminate between the two proposed 
models for IRS lA (jet versus equatorial wind). However, we note that the presence of the 
outflow in the 207200° direction (in addition to the outflow in the 1607340° direction), 
requires at least two independent powering YSOs in the region. In this way, within the 
equatorial wind scenario IRS lA would be the exciting source of the 160°/340° outflow (as 
proposed by Hoare 2006), while the new detected source IRS IC would power the 20°/200° 
outflow. 

Additional observations are required to discriminate in a more definitive way between 
the equatorial wind and the thermal jet scenarios. Any detection of large proper motions 
(hundred of km s^^) in the IRS ID and IE condensations in the northeast-southwest direction 
would favor the jet hypothesis. This test would require sensitive, high angular resolution 
continuum observations over several years. On the other hand, the study of the high velocity 
molecular gas in the vicinity of the source IRS lA could reveal the wide angle outflow 
expected from an equatorial wind. Any detection, with high angular resolution observations, 
of molecular gas correlated with the radio continuum source could trace the neutral part of 
the disk and also provide strong evidence favoring the equatorial wind scenario. In addition, 
as we have mentioned before, proper motion studies of the water masers observed in this 
region could be a powerful tool to discriminate between both models. On the other hand, 
with the present theoretical models (e.g. Drew ct al. 1998; Lugo et al. 2004), we cannot 
also discriminate between jet and equatorial wind models. However, we think that in an 
equatorial wind model the angular diameter might be less sensitive to the frequency of 
observation than in an jet model. In any case, more detailed theoretical studies will be 
necessary to address this important issue. 



5. Conclusions 

We have presented results of 1.3 cm continuum and water maser emission observations 
made with the VLA in its A configuration toward the SI 40 IRS region. We have also re- 
analyzed continuum observations at 0.7, 2, 3.5, and 6 cm. 

We observed IRS lA at all wavelengths (0.7, 1.3, 2, 3.5, and 6 cm) and also detected 
four new continuum peaks (IRS IB, IC, ID and IE). IRS IB is only detected at 1.3 cm, while 
IRS IC, ID and IE are detected at 3.5 cm. IRS IC, located- 0'/6 to the northeast of IRS lA, 
seems to be an independent source, while IRS ID and IE, located symmetrically with respect 
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to IRS lA, seem to be condensations associated with IRS lA. Under this scenario, IRS lA 
is not a single source, but a possible binary formed by IRS lA and IC. These two YSOs 
could explain the excitation of the two large bipolar molecular outflows observed toward 
S140 IRSl. We have also detected three water masers toward IRS 1 that most probably 
are tracing unbound motions and are associated with the bipolar molecular outflow in the 
northeast-southwest direction. 

In order to understand the nature of IRS lA, we have analyzed two scenarios: the 
thermal jet model and the equatorial wind model. The elongated morphology and spectral 
indices, a and j3, of the continuum emission in the 0.7 - 6 cm wavelength range, are consistent 
with IRS lA being a thermal jet. However, a photoevaporated disk (Hoare 2006, Gibb & 
Hoare 2007) could produce a similar morphology and spectral energy distribution. Proper 
motion measurements of the detected continuum clumps IRS ID and IE would help to 
discriminate between these two scenarios (jet versus equatorial wind). In addition, proper 
motion measurements of the detected water masers would also be very valuable in this sense. 

M.A.T. acknowledges the support from CONACyT grant 46157-E. J.M.T. acknowledges 
partial financial support from the Spanish grant AYA2005-08523-C03 
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Fig. 1. — Radio continuum contour maps showing consistently an elongated structure 
for S140 IRS lA at several wavelengths. Contours are -3,3,5,7,9,12,15 and 20 {panel a: 
1.3 cm; beam Cf.'ll x (y.'09), -3,3,5,7,9,11, and 13 {panel b: 2 cm; beam (y.'12 x (/.'lO), - 
3,3,5,7,9,12,15,20,30 and 40 {panel c: 3.5 cm; beam 0'/36x0'/31), and -3,3,5,7,9,12,15,20,30,40 
and 50 {panel d: 6 cm; beam 0'.'37 x 0'/31) times 200, 170, 43 and 40 //Jy beam"^, respec- 
tively, the rms noise of the maps. In all maps the beam is shown in the lower left corner. In 
order to make all continuum peaks of IRS lA coincide, we have applied an offset of (a, 5) 
= (-0f0017, -0'.'007), (+0f0017, +^.'009), and (+0f0016, +0'.'014) to the position of IRS lA 
at 2, 3.5 and 6 cm, respectively. The crosses show the position of the water masers detected 
in the region with their LSR velocity values (km s^^) indicated in panel b. The dashed and 
dotted lines in panel c indicate, respectively, the axes of the H2 jet and CO outflow seen on 
large scales. 
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Fig. 2. — Multiple radio continuum sources detected toward S140 IRS lA at different wave- 
lengths, a) 1.3 cm continuum contour map. Contours are -3,3,4,5,7,9,11,13 and 15 times 
230 //Jy beam~^, the rms noise of the map. The beam size is (X'OQ x (y/07. h) 3.5 cm con- 
tinuum contour map. Contours are -3,3,5,7,9,11,15,20,25 and 30 times 55 /xJy beam~^, the 
rms noise of the map. The beam size is 0'.'29 x 0'.'20. The crosses show the position of the 
water masers detected in the region. 
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Fig. 3. — Contour map of S140 IRS lA at 7 mm. Contours are -3,3,4,5,6 and 7 times 
310 //Jy beam~^, the rms noise of the map. The beam size (0'/04 x Qf!03) is shown in the 
lower left corner. The crosses show the position of the water masers observed toward the 
IRS lA. 



Table 1. Summary of the continuum observations toward S140 IRS 1 



Wavelength Observation HPBW Reference 
(cm) Date (arcsec) 



0.7 


1996 Nov 01 


~ 0.04 


VLA data archive (project AH5980) 


1.3 


1999 Jun 29 


~ 0.1 


This paper 


2 


1987 Sep 12 


~ 0.1 


Schwartz (1989) 


3.5 


1992 Nov 24 


~ 0.3 


Tofani et al. (1995) 


6 


1987 Sep 12 


- 0.3 


Schwartz (1989) 



-16- 




Fig. 4. — Flux density (top panel) and deconvolved angular size of the major axis (bottom) 
of IRS lA are plotted as a function of frequency. The solid lines represent the least-squares 
fits from which the spectral indices a and (3 are obtained for IRS lA. The filled triangles 
indicate the fiux density (upper panel) and the size of the major axis (lower panel) of IRS 1 
as obtained by Hoare (2006). These values were not used for the least-squares fits. The fiux 
density, S^, is given in mJy, the frequency, is given in GHz, and the deconvolved size of 
the major axis, Ojnaji is given in arcsec. 
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Table 2. Parameters of the continuum sources in S140 IRS 1 



Source Wavelength Position ^ Flux Density Angular Diameter 

(cm) Qi(1950) 5(1950) (mJy) 



IRS 


lA 


0.7 


22 


17 


41 


.083 


63 


03 


41, 


.65 


10.4±1.4 


0'.'09 X 0'.'02; P.A. 


61° 






1.3 


22 


17 


41 


.083 


63 


03 


41, 


.64 


9.1±0.7 


0'/18 X 0'/08; P.A. 


42° 






2 


22 


17 


41 


.085 


63 


03 


41, 


.67 


6.3±0.5 


0'.'27 X 0'.'03; P.A. 


45° 






3.5 


22 


17 


41 


.091 


63 


03 


41, 


.70 


3.9±0.1'' 


0'.'62 X O'.'ll; P.A. 


46° 






6 


22 


17 


41 


.087 


63 


03 


41, 


.69 


4.5±0.1 


0'.'68 X 0'.'07; P.A. 


41° 


IRS 


IB 


1.3 


22 


17 


41 


.070 


63 


03 


41, 


.57 


0.92^ 


< 0'.'08 




IRS 


IC 


2 


22 


17 


41, 


.136 


63 


03 


42, 


.08 


0.40^= 


< 0'.'15 








3.5 


22 


17 


41 


.128 


63 


03 


42, 


.24 


0.45^ 


< 0'.'2 




IRS 


ID 


3.5 


22 


17 


41 


.129 


63 


03 


41, 


.96 


0.38^ 


< 0'/2 




IRS 


IE 


3.5 


22 


17 


41 


.045 


63 


03 


41, 


.53 


0.38^ 


< 0'/2 





^Units of right ascension are hours, minutes, and seconds, and units of declination are 
degrees, arcminutes, and arcseconds. 

''This value is lower (~ 50%) than that reported by Tofani et al. (1995). 

'^Peak flux density. The beam size is 0'.'09 x 0'.'07 at 1.3 cm, 0'.'16 x 0'.'14 at 2 cm, and 
0'.'28 X 0'.'19 at 3.5 cm. 



Table 3. Parameters of the water masers observed toward S140 IRS 1 



Position ^ Vlsr Flux Density 

q;(1950) 5(1950) (km s-^) (Jy) 



22 17 41.099 63 03 41.64 7.7 0.35 
22 17 41.024 63 03 41.48 -16.0 2.29 
22 17 40.999 63 03 41.42 -13.4 1.44 



^Units of right ascension are hours, minutes, and sec- 
onds, and units of declination are degrees, arcminutes, 
and arcseconds. 



